Tunable Strongly Correlated Band Insulator 
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We introduce the notion of the strongly correlated band insulator (SCI) , where the lowest energy 
excitations are collective modes (excitons) rather than the single particles. We construct controllable 
1/N expansion for SCI to describe their observables properties. A remarkable example of the SCI 
is bilayer graphene which is shown to be tunable between the SCI and usual weak coupling regime. 
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Introduction - Idealized models of band insulators [1] 
are based upon determining the spectrum of the single 
particle excitations (SPEs), Fig. la (I), which are charac- 
terized by the fermionic statistics, momentum k, charge 
±e, spin 1/2. All the other excitations, i.e., electron-hole 
pairs are combined from the SPEs forming the particle- 
hole continuum [shaded region on Fig. la (II)]. 

Weak electron-electron repulsion does not change the 
gapped spectrum of the SPEs significantly, however it 
opens the decay channel of the SPE into three-particle 
continuum. This process is allowed only for the parti- 
cles with the energy e above the threshold, see Fig. lb 
(I) . The qualitative difference appears in the two-particle 
spectrum: electron-hole bound states (excitons) [2] are 
split down from the particle- hole continuum Fig. lb (II). 
Those discrete branches A'„(fc) (number of branches is in- 
finite for the interaction potentials lim^-^oo U{r)r^ = oo) 
can not decay unless their energies exceed some thresh- 
old. We will refer to this situation as weak-coupling insu- 
lator. All the thermodynamic and transport properties 
of such insulator are described by the SPE whereas exci- 
tons are responsible for the fine structure of the optical 
spectra. 

With the increase of the interaction the excitation hi- 
erarchy in the band insulator changes qualitatively, see 
Fig. 1 c. In this case, u > Xq{Q), and all the low tem- 
perature thermodynamics and the energy transport is 
contributed mostly by the excitons whereas the charge 
transport is determined by the SPEs. This leads to 
the different temperature dependence for the electric and 
thermal conductivities. Moreover, unlike in the case of 
the weak coupling insulators, the SPEs have a very nar- 
row stability range (Fig. Ic (I)), above which the elec- 
tron starts producing excitons similarly to well-known 
Schwinger mechanism [3] of vacuum polarization. For 
the same reason, only a finite number of the exciton 
branches are stable - all the other can decay into the 
two-exciton continuum (Fig. Ic (II)). We will call such 
system a strongly correlated insulator (SCI). 

In this Letter, we present the case study of the SCI 
using the bilayer-graphene (BLG) in a transverse electric 
field [4-6] as an example. A strong motivation for study- 
ing BLG in this context is the possibility to tune it from 



the weak coupling to the SCI, as it is discussed below. 
Using the number of electron species, A^ = 4, (two-fold 
valley and spin degeneracies) as a large parameter we ob- 
tain analytic results for observables determined by one- 
and two- particle excitations. We emphasize that our 
study is general and applicable to any SCI. 




FIG. 1. One- (I) and two- particle spectra (II) in the insulator, 
a) Non-interacting particles; b) the weak coupling; c) SCI; 
Particle-hole symmetry is implied for simplicity. 

Model- The band structure of BLG is described by the 
effective single-particle Hamiltonian [7] (/i = 1) 
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where uq is the interlayer asymmetry tunable by the 
transverse electric field, m is the effective mass, and 
P = (PxiPy) is the the Bloch momentum counted from 
the K point of the Brillouin zone. Excitations whose mo- 
menta are near K' are described by the parity conjugate 



of the Hamiltonian (1) and form another particle species 
with identical properties. The four-fold degenerate spec- 
trum of the SPEs corresponding to Eq. (1) is 

eik) = {ul + fcV(4m2))'/' « uo + ky{8u„m^). (2) 

We will see that SCI behavior of BLG is associated 
with the long range part of the Coulomb interaction 
v(r) = e^/r so that the effective Hamiltonian is 

H = (frt/j'^Ho^ + - d^rcfr' p{y)v{y - v')p{y'). (3) 

We omitted here the short range interaction with non- 
trivial matrix structure. Those terms are related to the 
symmetry breaking in BLG and get renormalized with 
increasing the linear scale [8]. All those renormalizations, 
as well as the renormalisation of the electron mass and 
interlayer asymmetry, stop at distances of the order of 



(a) 



(. = l/^?71Uo, 



(4) 



and are incorporated into the SPE spectrum (1). In 
Eq. (3), Hq given in Eq. (1) operates on a two-component 
vector ■0(r) = [-0yi(r),^B(r)]'^ in the sublattice space 
(summation over the spin and valley indices is implied), 
p =: ^T^ • jg lY^Q normal-ordered particle density. 

Dimensional analysis of the Hamiltonian (3) reveals 
only one dimensionless coupling constant 



K ~ me /uq, 



(5) 



which diverges a.s uq —>■ signaling the SCI. 

Screened interaction in large N approximation is ob- 
tained as a resummation of Fig. 2 a: 



V{q, iio) ^ v{q) [I + Nv{q)n{q, ito)] 



(6) 



where v{q) = 27re^/g and Il{q,iuj) is the polariza- 
tion function in the momentum-frequency domain with 
asymptotic expressions shown in Fig. 2 b). Taking iV ^ 1 
limit, one finds k to drop out of the expression giving 



V{q,iuj)^ l/[NIl{q,iLu) 



(7) 



Thus, the perturbation series becomes an expansion in 
powers of 1/A^. Such an expansion, however, contains 
singular terms caused by H(g — > 0, iuj) — > 0. The expres- 
sion for the screened static potential is most instructive: 



V{r) 



d^qe' 



-y(q,o) 



2^1n^, 

A' „ r ' 



(27r)2 

where new spatial scale appeared: 
i?, = kNE. 



r <^ R^ 
r:S>R* 



, (8) 



(9) 



Notice that V{£) > ua at at small enough uq, which 
makes the SCI insulator regime possible. At larger uq 
interaction is weak and may be treated pcrturbatively. 
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FIG. 2. a) Resummation of the leading A^ loops leads to the 
RPA screening of the interaction potential, b) Asymptotic 
regions for the polarization operator, c) The leading-order self 
energy correction, d) Resummation of the leading logartihm 
using Ward-Takahashi identity (e). 



SPE spectrum and particle residue - We start with 
calculating the leading in 1/A^ self-energy correction. 
Fig. 2c. Finding the pole of the resulting Green func- 
tion, we obtain instead of Eq. (2) 



eik) 



0.680 k^ 



N 2m Sm'^uo 
The single particle gap is strongly enhanced: 



il + 0{l/N)). (10) 
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(11) 



In the limit of k — > oo this enhancement diverges which, 
as we will see shortly, signals the transition to the SCI. 
The physical meaning of gap renormalization is that the 
band insulator can not completely readjust itself when a 
charged particle is introduced. The interaction of the ex- 
tra particle with the dipoles from the distances between i 
and i?* leads to the logarithmic divergence for i?* — > oo. 
The second term in Eq. (10) is the interaction-induced 
negative curvature analytic in 1/A'^. It significantly ex- 
ceeds contribution of Ref. [7] to the resulting "mexican 
hat" spectrum shown on Fig. Ic (I). 

The single-particle residue found from the same dia- 
gram of Fig. 2c is given by 



Z Kl- (S/tt^N) In(KiV) In [A/uq) , 



(12) 



where A w 0.2 eV is the upper limit for the applicability 
of the two-band model (1). 

Equations (11) and (12) contain logarithmically diver- 
gent factors, which makes it necessary to sum up all or- 



dcrs of the perturbation theory in 

A= {lnKN)/{27rN). 



(13) 



The summation procedure outhned below shows that the 
gap (11) does not acquire any higher order in A correc- 
tions whereas Eq. (12) changes to 

Z = exp [- (S/TT^iV) \n{nN) hi (A/mq)] ■ (14) 

Let us sketch derivations of Eqs. (11) and (14). We 
notice that the dominant logarithmic divergence is con- 
tributed by the momenta transfer q^ ^ 1 through any 
interaction propagator. This allows for following simpli- 
fications: (i) to factorize the polarization operator as 



U{q, iuj) = q'^/ [uag{iuj/uo)] , 



(15) 



where g{z) is a function well defined in the limit n — > oo; 
(ii) to neglect the change in the momentum k of any 
single particle Green function (GF). The latter allows 
the integration over the momenta in each interaction line 
separately: the screened interaction (6) is replaced with 



e<i 



V{q, iuj)(fiq/ {2Tr)'^ = uoXg {iuj/uq) 



(16) 



A non-trivial aspect of Eq. (15) is the appearance on the 
bare scale uq rather than the rcnormalized gap u in the 
right-hand-side of the equation. The reason for this is 
that divergence ex InA'^K can not enter into II, because 
the screening is determined by the neutral dipoles [9]. 

Dyson equation for the resulting strong coupling the- 
ory, Fig. 2d, is then closed using the Ward-Takahashi 
identity Fig. 2e. This results in exact time ordered GF 



/oo 
dte'''Go{Kt)e-^^^'\ (17) 
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pA/uo 

F{t) = ig{Q)u\t\ + dzh- e-*"l*l j W{z). 



where Gg (k, t) is the bare GF in the momentum-time 
domain, and W{z) = i [g{z + iO) — g{z — iO)] /{nz^). 

Expression (17) is valid for any function g. For the po- 
larization loop neglecting vertex corrections and renor- 
malization of GF, N -^ oo [9], we find g{0) = 6tt and 



W{z) 



IGttz 



In 



z + 2 



4z 



z^ +4 



Then, F{t) = 67rm|i|+(16/7r) ln(A/uo) and Eq. (17) leads 
to Eqs. (11) and (14). 

Equations (17) enables one to calculate not only the 
particle pole but also the incoherent contribution describ- 
ing the coupling of the extra electron introduced or ex- 
tracted from the system (as in the tunneling or photoe- 
mission experiments) with many particle continuum. To 



A{e) 




FIG. 3. The electron spectral weight A{e) at different values of 
the coupling constant At. For clarity, the curves corresponding 
to ten equidistant values of k = 4, 8, 12, . . . , 40 are vertically 
offset and the quasiparticle peak is artificially broadened and 
scaled down by a factor of 100. 



do so, we evaluate A{e) = (l/7r)|ImTrG(0, e)| by numer- 
ical integration of Eq. (17). (Due to the electron-hole 
symmetry, ^(e) = A{—e), and the structure and finite k 
is similar). The result plotted in Fig. 3 shows the sin- 
gle particle peak at e = u ^ uq and the threshold at 
e = u + 2uo due to the coupling to the three-particle 
continuum [10] 

Exciton spectra - Next we discuss the collective spec- 
trum of the system, revealed as the poles of two-particle 
propagators of the system. To calculate those poles, X, 
in leading 1/A^ approximation it is sufficient to neglect 
the retardation in the interaction potential (6) and con- 
sider Schrodinger's equation for an electron and a hole. 



Pe 



-v\ 



+ {2u - V(r)\ \ ^(r,,rn) - X^{t,, m). (18) 



with V(r) defined in Eq. (8), and r = rg — v^. Unbound 
states with X > 2u correspond to the particle-hole con- 
tinuum and the bound states are the exciton lines. 

Due to the non-parabolicity of the one-particle spec- 
trum the motion of the exciton center of mass P = 
Pe + P/t can not be separated from the relative motion, 
and we consider here the case of P = relevant for optics. 
In this case, the levels are labeled by the normal angu- 
lar momenta |j| = 0, 1, 2, . . . and by the radial quantum 
number n = 0, 1, ... . Dimensional analysis gives the low 
lying states with the size of the order of ^ ^ i?*, thus 
[2u - V(t)\ == uo [2 + (3/2Ar) Inr/^] does not depend on 
K ^ 1, and, for ?i + |j| < A^k, 



Xi ^Uq 



AT \4 A^ 






(19a) 



where ^^ arc the eigenvalues of the differential operators 

h? = {-d^ / dz^ - z-'^d/ dz + f / z^f +\n z. (19b) 

Results of the numerical diagonalization of h^ are plotted 
on Fig. 4. Note that the levels with X-j-^ > 2Xq decay into 
two exciton continuum as in Fig. Ic (II). 
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FIG. 4. The energy levels of the exciton in units of uq for 
N — 4. Hatched region is the two-exciton continuum. Side 
panels are the sketches of the absorption (left) and Raman 
spectrum (right). Bottom table are the oscillator streghts of 
the transitions in units of the strongest line, e^^o are the com- 
plex polarization vectors for the incoming (i) and outgoing 
(o) light, Q is the base line frequency and 71 ~ 0.2 eV. The 
coefficient 1/3 ~ {ShiTr/u)^ is not zero only due to the trigonal 
warping of the spectrum, where Slitt — 1 meV is the energy 
of the Lifshitz transition in gapless, uq — 0, BLG [7, 8]. 



Since the spectrum of the excitons (19a) does not de- 
pend on the coupling constant k, the system can be quite 
easily tuned in the SCI regime, u ^ Xq. 

The fact that energies of the lowest-lying excitons in 
Eq. (19a) scale with bare gap parameter uq rather than 
the renormalised single-particle gap u justifies the scal- 
ing form of the polarization operator Eq. (15). Indeed, 
for the strong coupling regime, the polarisation of BLG 
occurs via appearence of virtual excitons with j = ±1, 
so that n(g <^ ^~^,w) is generic for both weak coupling 
regime and SCI. Moreover, g{0) acquires only 1/iV cor- 
rections, so that Eq. (11) does not change. The fine struc- 
ture of g{z) shows the exciton resonaces it changes the 
threshold in ^(e) by the exciton binding energy ~ uq/N 
and introduces an additional fine structure which may be 
distinguished in higher derivatives of A{e) 

Exciton lines in SCI optics - The exciton lines (19a) 
are degenerate {N"^ for j=0, and 2iV^ for \j\ > 0) due 
to the iV-fold degeneracy of the electron spectrum. Such 
degeneracy is lifted due to the crystalline symmetry [11]. 

First, the 4-fold spin degeneracy is split to the S ~ 
singlet and S = 1 triplet states due to the exchange in- 
teraction. In the absence of spin-orbit interaction triplet 
excitons can not be observed in optical experiments. 
Spin-singlet states are further split due to the trigonal 



symmetry of the bilayer crystal and should be classified 
according to the irreducible representations of its pla- 
nar group, only Ai^2 and E representations are optically 
active. Such a classification is presented in Fig. 4 to- 
gether with the selection rules (^1 is not active in Ra- 
man because of the electron- hole symmetry) . These rules 
for the bright i5-exciton absorption (as well as lumines- 
cence), a^ -^ X^ — >• cr^ are determined by the form of 
the interband current operator derived from the Hamil- 
tonian Eq. (1), and a trigonal warping term [7] due to 
skew intcrlaycr hopping, for a weak transition a —>■ X^. 
The selection rules for Raman processes are determined 
by the electron-two photon interaction via virtual in- 
termediate state [12], with the dominant transition A2: 
a^ -^ ±i(j^ + X"^ and a satelitc E: a^ ^ a^ + X^ . 

In conclusion, we presented a general and control- 
lable theory of a strongly correlated insulator (SCI): a 
band insulator where the the spectrum of excitons lies 
deep below the lowest branch of the single-particle spec- 
trum. Gapped bilayer graphene is not the only example 
with such properties, the list of other potential SCIs in- 
cludes quantum wells of semimetal compounds, such as 
Bii_2;Sba;, or silicene [13] in a transverse electric field 
[14]. However, BLG is unique in its tunability from the 
weak coupling to the SCI. In the lab, this tuning can be 
achieved by the application of the electric field normal to 
the BLG plane. 

We acknowledge discussions with L.I. Glazman and 
support by US DOE contract No. DE- AC02-06CH11357 
(lA), EPSRC EP/G041954/1 (VF), ERG (VC), and the 
Royal Society (VF). 



[i; 
[2: 

[3] 
[4] 
[5; 

[6; 

[7 

[8 
[9; 

[lo; 

[11 

[12 
[13 
[14 



See e.g. C. Kittel, Introduction to Solid State Physics, 
8th edn, Wiley Press. 2005. 

G. Wannier, Phys. Rev. 52, 191 (1937); N.F. Mott, 
Trans. Faraday Soc. 34, 500 (1938). 
J.S. Schwinger, Phys. Rev. 82, 664 (1951). 
K.S. Novoselov, et ai, Nature Phys. 2, 177 (2006). 
Y. Zhang, et ai, Nature 459, 820 (2009). 
R.T. Weitz, et ai, Science 330, 812 (2010). 
E. McCann and V.I. Fal'ko, Phys. Rev. Lett. 96, 086805 
(2006). 

Y. Lemonik, et al, Phys. Rev. B 82, 201408 (2010). 
We will elaborate on this after the exciton is considered. 
In fact A{e) is singular at e„ = u + 2nuo but those addi- 
tional singularities are seen only in higher derivatives. 
Excitons: selected chapters, E.I. Rashba and M.D. Sturge 
eds, North-Holland Pub. Co. (1987). 
M. Mucha-Kruczynski, O. Kashuba, and V.I. Fal'ko, 
Phys. Rev. B 82, 045405 (2010). 

S. Cahangirov, et al, Phys. Rev. Lett. 102, 236804 
(2009). 

V. Zolyomi, N. Drummond, I. Aleiner, and V. Fal'ko - 
unpublished 



